Astronomy & Astrophysics manuscript no. merger-rate- VVDS-arxiv2 © ESO 2009 

February 6, 2009 



G\ 



Received 11/07/2008; accepted 31/01/2009 



The VIMOS VLT Deep Survey 

Evolution of the major merger rate since z ~ 1 from spectroscopically confirmed 

galaxy pairs. * 

L. de Ravel 1 , O. Le Fevre 1 , L. Tresse 1 , D. Bottini 2 , B. Garilli 2 , V. Le Bran 1 , D. Maccagni 2 , R. Scaramella 413 , M. 

Scodeggio 2 , G. Vettolani 4 , A. Zanichelli 4 , C. Adami 1 , S. Arnouts 23 ' 1 , S. Bardelli 3 , M. Bolzonella 3 , A. Cappi 3 , S. 

Chariot 810 , R Ciliegi 3 , T. Contini 7 , S. Foucaud 21 , R Franzetti 2 , 1. Gavignaud 12 , L. Guzzo 9 , O. Ilbert 20 , A. Iovino 9 , F. 

Lamareille 7 , HJ. McCracken 10,11 , B. Marano 6 , C. Marinoni 18 , A. Mazure 1 , B. Meneux 22 - 24 , R. Merighi 3 , S. 

<N Paltani 1516 , R. Pello 7 , A. Polio 117 , L. Pozzetti 3 , M. Radovich 5 , D. Vergani 2 , G. Zamorani 3 , E. Zucca 3 , M. Bondi 4 , A. 

^O ■ Bongiorno 22 , J. Brinchmann 19 , O. Cucciati 9 ' 14 , S. de la Torre 1 , L. Gregorini 4 , P. Memeo 2 , E. Perez-Montero 7 , Y. 

r <P \ Mellier 10 - 11 , P. Merluzzi 5 , and S. Temporin 9 

Mm . 

_ , (Affiliations can be found after the references) 

Y*\ ABSTRACT 

o 

5-H Context. The rate at which galaxies grow via successive mergers is a key element to understand the main phases of galaxy evolution. 
(y-j Aims. We measure the evolution of the fraction of galaxies in pairs and the merging rate since redshift z ~ 1 assuming a (H = 
C3 10kms~ l Mpc~' ,Q M = 0.3 and £2a = 0.7) cosmology. 
1 ' . Methods. From the VIMOS VLT Deep Survey we use a sample of 6464 galaxies with I AB < 24 to identify 314 pairs of galaxies, each member 

with a secure spectroscopic redshift, which are close in both projected separation and in velocity. 

^N Results. We estimate that at z ~ 0.9, 10.9 ± 3.2% of galaxies with M B (z) < -18 - Qz(Q = 1.11) are in pairs with separations Ar p < 20h~ l kpc, 

i-** Av < 500km/s, and with AM B < 1.5, significantly larger than 3.8+ 1.7% atz ~ 0.5; thus, the pair fraction evolves as (1 +z) m withm = 4.73 ±2.01. 

For bright galaxies with M B (z = 0) < -18.77, the pair fraction is higher and its evolution with redshift is flatter with m = 1.50 ± 0.76, a property 

also observed for galaxies with increasing stellar masses. Early-type pairs (dry mergers) increase their relative fraction from 3% at z ~ 0.9 to 12% 

at z ~ 0.5. The star formation rate traced by the rest-frame [Oil] EW increases by 26 ± 4% for pairs with the smallest separation r p < 20h~ [ kpc. 

Following the prescription to derive merger timescales of Kitzbichler & White (2008) we find that the merger rate of M B (z) < — 18 - Qz galaxies 

|^! , evolves as N mg = (4.96 ± 2.07) x 10~ 4 ) x (1 + z) 220±0J1 mergers Mpc^Gyr 1 . 

f*^ • Conclusions. The merger rate of galaxies with M B (z) < - 18 -Qz has significantly evolved since z ~ 1 and is strongly dependent on the luminosity 
ryQ ' or stellar mass of galaxies. The major merger rate increases more rapidly with redshift for galaxies with fainter luminosities or stellar mass, while 
f"~) the evolution of the merger rate for bright or massive galaxies is slower, indicating that the slow evolution reported for the brightest galaxies is not 
• • ■ universal. The merger rate is also strongly dependent on the spectral type of galaxies involved. Late-type mergers were more frequent in the past, 
^ ■ while early-type mergers are more frequent today, contributing to the rise in the local density of early-type galaxies. About 20% of the stellar mass 
kj ' in present day galaxies with Iog(M/M e ) > 9.5 has been accreted through major merging events since z = 1. This indicates that major mergers have 
r% ' contributed significantly to the growth in stellar mass density of bright galaxies over the last half of the life of the Universe. 

" Key words. Galaxies: evolution; Galaxies: interactions 

1. Introduction the two components have more or less the same mass, disks 

could be transformed into spheroidals, as predicted using de- 



In the current hierarchical structure formation paradigm, the tailed simulations dCombes F„ 2004; Mihos & Hernquist, 1996 



mass assembly in galaxies proceeds via a process of coa es- | Con selice et al., 2006| . It is also expected that major merger 

cence between increasingly more massive dark matter halos. ' eV ents profoundly modify the spectrophotometry properties of 

This halo merging tree history can be quantified by a halo ^ laxies involvedi for instance tr i gge ring a burst of star 

merger rate, measuring the growth of mass per average mass in formation ( |p a tton et al., 2 005). Galaxies in the process of 

a representative volume of the Universe. However, these mod- { are ob ' served , however me contribution of this pro- 

els do not directly predict a growth of galaxy mass via mergers cess tQ ^ evolution of the lobal j population is not yet 

(Moore etal.,200g, and the actual contribution of mergers to isd constrained Indirect evidence for m is also in _ 

the evolution of galaxies remains poorly predicted. ferred from Qther j parties like the luminosity or mass 

Merging two galaxies is potentially a very powerful process. function The luminosit of the red bul dominated popu i a . 

It is possible that during m ajor merger events, i.e. mergers where tkm of galaxies is measured to increas e since z ~ 1 p art 

„ , r,. . . . . T , ,, , , . , ,„ c of which could be produced by mergers (Ilbert et al., 20061. it 

Send offprint requests to: L. de Ravel : loic.deravel@oamp.lr r . J . ° \ : —4 ' 

* , , , ,. , . . , ... ., „ „ ., seems possible that the increase in the density of lntermedi- 

based on observations obtained with the European Southern r . . J 

Observatory Telescopes at the Paranal Observatory, under programs ate mass early-type galaxies since z ~ 1 may be happening 

072 A-0586 and 073 A-0647 at me ex P ense of late-type galaxies involved in merging events 
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( |Tresse et~ al., 2007). Merging is therefore potentially a very im- 
portant physical phenomenon which could drive the evolution 
of galaxies along cosmic time. The average numbers of merger 
events needed to build a typical M* galaxy, the contribution 
of mergers to the mass growth of galaxies, or the identifica- 
tion of a prefered time in the life of the Universe when merg- 
ers were more frequent, are all important elements to help to- 
wards our understanding of galaxy evolution. It is then crucial 
to quantify the contribution of merging to the evolution pro- 
cess and its impact on important quantities like the cosmic star 
formation rate (e.g. |Bouwens et al., 2006, |Tresse et al., 2 007, 
Woods et al., 2006]) or the global stellar mass density (e.g. 
Arnouts et al„ 2007HBundy et al, 2005J|Pozzetti et al„ 20071 ). 



To estimate the contribution of mergers to the formation and 
evolution of galaxies is not a trivial task. In the nearby Universe 
merger events can be identified aposteriori from perturbed mor- 
phologies, wisps, tails, and other peculiar signatures seen at low 
surface brightness. Only recently volume complete measure- 
ments of the merger rate in the nearby Universe are becoming 
available. In the Millennium catalogue, especially tailored to a 
volume complete identification of merging events, de Propris 
et al. (2007) use the relative velocity measured from spectro- 
scopic redshifts to confirm true galaxy pairs in the process of 
merging. They find that the merger fraction is 2% at a mean 
redshift of 0.06, refining earlier estimates based on pair frac- 
tion ( [Patton et al., 20001 [Patton et al., 2002) ). At higher redshifts, 
searching for evidence for past mergers becomes increasingly 
difficult, because the residual signatures of mergers often have a 
too low surface brightness. At redshifts z > 0.3, it is therefore 
easier to search for 'apriori mergers', encounters that are likely 
to lead to a merger event, rather than to look for 'aposteriori' 
signs of past mergers. When two galaxies are close together in 
space, and depending on their relative velocities, gravity is acting 
to bring them closer for a bound system that will merge. A mea- 
sure of the merging frequency is then to count galaxy pairs with 
a separation and velocity difference such as they are likely to be 
gravitationally bound and destined to merge. By selecting pairs 
of galaxies with similar magnitudes and hence approximately 
with similar masses, one can focus on major merger events. 
They are able to significantly contribute to the mass assembly, 
to modify morphologies, as well as to significantly alter the star 
and gas content of the incoming galaxies. Assuming that a dy- 
namically bound system of two galaxies will most likely evolve 
into one more massive galaxy, one can then derive the merger 
rate from the pair count. A major uncertainty of this estimator is 
the timescale upon which a merger will be completed. N-body 
simulations are then used to provide reasonable estimates of the 
merger timescales (Conselic e et al., 20061 Kitzbichler & White, 
2008). 

Cold dark matter simulations show that the evolution 
of the dark matter halos merger rate follows a power law 
N mg = N mgt o(l +z) m where N mg fl = N mg (z - 0) is the local value, 
and m parameterizes the evolution. While some simulations 
predict that m should have 2.5 < m < 3.5 (Gottl0ber, 2001), 
measuring m directly from galaxy samples is an important step 
to understand the evolution of galaxies. Many observational 
attempts have been carried out to track the evolution of the 



merger rate as a function of redshift (e.g. Burkey et al., 1994 
Carlberg, Pritchet & Infante, 1994[ |Yee & Ellingson, 1994 

IPatton et a~2 000 



Patton et al., 1997 



Patton et al., 2002 



|Le Fevre et al, 2000 



Conselice et al., 2003, Lin et al., 2004 



Kartaltepe et a l., 2007[ |Lin et al„ 20081 |Lotz et al., 2 008 



Kampczyk et al., 2007 1. Even though, m remains poorly con- 
strained with < m < 6, meaning either no evolution of the 



merger rate with cosmic time, or a strong evolution. Part of 
this large range of values can be understood as coming from 
the different criteria used to identify merger candidates, or the 
photometric band used to identify pairs (Bundy et al., 2005). 
Furthermore, comparing measurements at low and high redshifts 
from different surveys is complicated due to the different selec- 
tion functions used. At redshifts z > 0.3, most pair counts so far 
have been performed from a measurement of the number of pairs 
observed on deep images, with either a photometric redshift 
of the galaxies (e.g. |Conselice et al., 2003) |Lotz et aLT"2 008), 
or a spectroscopic redshift of one of the galaxies in the pair 
(e.g. Patton et al., l"997l |Le Fevre et al., 20001 ). The effect of 
contamination by galaxies projected along the line of sight 
producing false pairs is then estimated from galaxy counts, and 
the observed pair fraction is corrected to get an estimate of 
the true pair fraction. As redshift increases, projection effects 
become increasingly important making it difficult to estimate 
the true pair fraction, creating a fondamental uncertainty in the 
measurement of m. At z ~ la galaxy with a luminosity L* has 
a 40% probability to have a galaxy with a similar magnitude 
but at a different redshift projected within an apparent radius of 
2Qh~ l kpc flLe Fevre et al., 20001 ). 

To overcome these limitations, the most secure method to 
identify a physical pair of galaxies is to obtain a velocity 
measurement of each galaxy in the pair, enabling to identify 
pairs of galaxies which are most likely to be gravitationnaly 
bound. Only recently samples with spectroscopic redshifts for 
both galaxies in a pair are becoming available ( ]Lin et al~2 007 
|Lin et al., 20081 ). I n this paper we use for the first time a com- 
plete redshift survey to z ~ 1 and as faint as Iab - 24 to 
securely identify pairs with both galaxies having a spectro- 
scopic redshift. We use the VIMOS VLT Deep Survey (VVDS) 
( [Le Fevre et al., 2005al ), to search for galaxy pairs and to derive 
the pair fraction and the merger rate evolution. We present the 
galaxy sample and the methodology to build a pair sample in 
Section 2, we derive the pair fraction evolution in Section 3, 
and we examine the spectrophotometric properties of galaxies 
in pairs in Section 4. We compute the merger rate in Section 5. 
We evaluate the fraction of the stellar mass involved in mergers 
since z ~ 1 in Section 6, and conclude in Section 7. We adopt 
a Ho = 70kms~ l Mpc~ l , Q.a = 0.7 and Q„, = 0.3 cosmology 
throughout this work and magnitudes are given in the AB sys- 
tem. 



2. Identification of galaxy pairs 

2.1. VVDS overview 

We use the deep sample from the VIMOS VLT Deep Survey 
on the 0216-04 field. Data have been obtained with the Visible 
Multi Object Spectro graph (VIMOS) on the ESO-VLT UT3 
flLe Fevre et al., 20031 ). A total of 9842 objects have been ob- 
served in the VVDS-Deep field over a total area of ~ 0.5 deg 2 , 
selected solely on the basis of apparent magnitude 17.5 < Iab ^ 
24. The mean redshift of the sample is z - 0.76. The veloc- 
ity measurement of each galaxy redshift has an accuracy of 
~ 276km/ s (Le Fevre, 2005a). A strategy of multiple spectro- 
graph passes have been used (Bottini et al., 2005), defining ar- 
eas where targets have been randomly selected in four separate 
observations, and another area where two independent observa- 
tions have been performed, leading to an effective random sam- 
pling of the galaxy population of respectively ~ 35% and ~ 20% 
for each of these two areas (with respectively Sn p = 0.17 deg 2 
and Sip = 0.32 deg 2 ). We use a catalogue which contains 6464 
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objects in an effective area of ~ 0.5 deg 2 , using only the most 
secure redshifts, i.e. quality flags 2,3,4 and 9 for primary targets 
and 22, 23, 24 and 29 for secondary targets. Flags 2 ,3 ,4 cor- 
respond to redshifts measured with a confidence level of 80%, 
95% and 100%, respectively, and flag 9 indicates spectra with 
single emission line (see |Le Fevre et al" , 2005a for details). 

2.2. Selection of pairs in the VVDS 

We have identified pairs in two ways. First, we have searched 
in the main VVDS catalogue to find pairs of galaxies close in 
separation perpendicular to the plane of the sky using the an- 
gular distance at the redshift of the pair, and close in veloc- 
ity along the line of sight as derived from the redshift mea- 
surements. Secondly, we have visually examined the 2D spec- 
tra to identify secondary objects close to a primary VVDS tar- 
get which have escaped the automated spectra detection algo- 
rithm (Scodeggio et al., 2005) because their angular proximity 
to the main target along the slit creates a blend of the two spec- 
tra at the faint isophotes used for detection. We looked for ev- 
idence for 2 continuum traces next to each other, with a clear 
separation of the objects profile along the slit. The ID spec- 
trum of the companion was then extracted and its redshift mea- 
sured using the cross-correlation with templates as done for 
the main VVDS sample , and was assigned a flag 3X, with 
X following the flag nomenclature of the survey as described 
in Section 2.1. We then search in the parent VVDS imaging 
catalogue for the object responsible for the secondary trace, 
providing its sky coordinates, the magnitudes and colors. In 
case the photometric catalogue did not identify the companion 
also because of blending, ugri and z images from the CFHTLS 
(http://www.cfht.hawaii.edu/Science/CFHLS/) were examined 



and the multi-band photometry of the companion was performed 
using flux extraction in image areas isolating the object. This 
process concerned mainly objects with separations 1 < 9 < 2 
arcseconds. 

Our final catalogue contains all primary target galaxies with 
secure redshift measurements (flags 2, 3, 4, 9) and Iab < 24 
(6287 objects), all secondary target galaxies identified by the au- 
tomated spectra extraction program (flags 22, 23, 24, 29) (160 
objects) and all the companions identified through our visual ex- 
amination of 2D spectra (flags 32, 33, 34, 39) (17 objects). 

To create the pair catalogue, we first compute two quantities: 
the projected separation r p and the line-of-sight velocity differ- 
ence Av. For a pair of galaxies with redshift z, and Zj and an 
angular separation 6 these parameters are given by : 



6 x d A (z,„), where z„ 



Zi + Zj 



(1) 



Av - c 



\Zi 



i+z n 



where dA(z m ) is the angular diameter distance at the mean pair 
redshift z m , and c the speed of light. 

The I AB < 24 selection of the VVDS implies that galaxies in 
the sample have an absolute magnitude Mb < -19.11 at z = 1. 
We are missing pairs for which one member of the pair is fainter 
than this limit, for which we need to apply a completeness cor- 
rection as described in Section 2.4. From the luminosity function 
of the complete VVDS sample, we know that the characteristic 
luminosity M* in B-band evolves with redshift as Q(z) - 1.11 Xz 
filbert et al., 2005[ ). We have therefore applied a magnitude evo- 
lution Mb — -18 - Q(z) to our absolute magnitude cutoff when 
looking for pairs (see Figure[T). 



^3 



a 

C3 



^3 

< 




0.4 0.6 

redshift 

Fig. 1. Absolute B band magnitude distribution versus redshift 
of the primary galaxies in a pair with r p < 100/r'fcpc and 
Av < 500km/ s (filled symbols) compared to the underlying pop- 
ulation of galaxies (open symbols). The line indicates the limit in 
absolute magnitude used to identify pairs with M B < -18 — Q(z). 



2.3. The VVDS pair catalogue 

Using equation (1), we have identified 702 simple pairs and 190 
triplets with r" p wx = ISOh^kpc and Av max = 2000 km/s. To 
select major mergers, we have imposed an absolute magnitude 
difference between the two members of a pair in the B band of 
AMb < 1.5 mag (see Figure |2). On the left panel of Figure [3] 
we present the number of pairs with Av mc,x - 2000 km/s and 
AMb ^ 1 .5 as a function of r p and on the right panel, the number 
of pairs with r'™ x - 150h~ l kpc and AM B < 1.5 as a function 
of Av. 

The two-point correlation function £(r) describes the ex- 
cess probability of finding a galaxy at distance r from a galaxy 
selected at random over that expected in a uniform, random 
distribution. This function is usually parametrised by a power 
law with correlation length r : £(r) = (r/ro)~ y . Integration 
over this function yields a number of pairs that varies as r p Y 
(Patton et al, 2002). Using the mean slope y ~ 1.7 of the corre- 
lation function found in the VVDS ( |Le Fevre et a l, 2005b), we 
expect an increase of the number of pairs ~ r p :i . This is in good 
agreement with our pair counts which gives a slope of ~ 1 .24 as 
shown in the left panel of Figure [3] 

We identify the number of pairs as a function of separations 
r max and Av max in Ta ble[T] For r p < 20/T 1 kpc, Av < 500 km/s 
and imposing at least one of the pair members to have M B < 
-18-<2(z), we have a total of 36 pairs. The fraction of close pairs 
added by the visual examination of the 2D spectra is ~ 10% for 
pairs with a projected separation less than 2 arcseconds. 



Table 1. Number of pairs with AM B mx = 1.5 and M B < -18 
Q(z). In brackets the number of triplets is given. 

20/T 1 kpc 30k- 1 kpc 50/r 1 kpc 100ft-' kpc 



500 km/s 
1000 km/s 
2000 km/s 


36(0) 
48(0) 
50(0) 


51(0) 
68(0) 
73(0) 


102 (6) 
131 (10) 
143 (13) 


202 (22) 
267 (32) 
314(46) 



We give the list of all 36 pairs with r p < 20h ' kpc, Av < 
500 km/s, AM B < 1.5 and M B < -18 - Q(z) in Tablefjand we 
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Fig. 2. The line-of-sight velocity difference Av as a function the projected separation r p for all pairs with r p < 1 50 h l kpc and 
AM" MV = 1.5 mag. 
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Fig. 3. Left : Number of pairs as a function of r™ flJC for all selected pairs with Av < 2000 km/s and AMb < 1.5 with M# < 



-18 - Q(z). The line represents the best fit to the number of pairs with N pa i rs oc r , comparable to the expectation oc rb 3 from the 
angular two-point correlation function. Right : Number of pairs as a function of Av max for all selected pairs with r p < 150/z L kpc 
and AM B < 1.5 with M B < -18 - Q(z). 



display the postage stamps and spectra of each of these pairs in 
Figure HI 
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Table 2. List of the 36 spectroscopic pairs with rJ" 0A = 20h~ l kpc, Av max - 500 km/s and AM™ ax = 1.5 mag selected in the bright 
Mb(z - 0) < -18 sample. Id's beginning by "p" are manually extracted. Pair numbers can be use to retrieve postage stamps and 
spectra in Figure|4] R. A. (2000) and Dec. (2000) are given for the first galaxy. 



Pair number 


Idl 


Id2 


R. A. (2000) 


Dec.(2000) 


zl 


z2 


Z-mean 


r p (h { kpc) 


Av (km/s) 


AM B 


en 


1 


020236244 


020236318 


36.533573 


-4.483242 


0.9252 


0.9266 


0.9259 


12.0 


217.9 


0.16 


2.2 


2 


020260617 


020260636 


36.720883 


-4.429412 


0.2657 


0.2658 


0.2657 


12.8 


23.7 


0.02 


4.5 


3 


020260559 


020260085 


36.655374 


-4.427911 


0.7088 


0.7085 


0.7087 


17.4 


52.6 


1.44 


3.5 


4 


020274782 


020273998 


36.521896 


-4.396756 


0.6295 


0.6320 


0.6308 


20.0 


459.6 


1.04 


4.2 


5 


020314240 


020314107 


36.585748 


-4.309167 


0.6888 


0.6872 


0.6880 


10.4 


284.2 


0.92 


2.1 


6 


020461143 


020461037 


36.699123 


-4.399960 


0.7047 


0.7043 


0.7045 


15.9 


70.4 


0.22 


3.2 


7 


020199214 


020199508 


36.277674 


-4.557759 


0.9126 


0.9107 


0.9116 


15.0 


298.0 


0.59 


2.7 


8 


020205594 


020204675 


36.714690 


-4.545729 


0.6309 


0.6323 


0.6316 


13.6 


257.2 


0.11 


2.8 


9 


020208200 


020207985 


36.647773 


-4.538258 


0.6946 


0.6929 


0.6937 


13.6 


300.9 


0.98 


2.7 


10 


020231154 


020230801 


36.518347 


-4.493709 


0.9265 


0.9250 


0.9257 


8.0 


233.5 


0.73 


1.5 


11 


020234145 


020234032 


36.685092 


-4.486811 


0.8859 


0.8885 


0.8872 


17.4 


413.0 


0.95 


3.2 


12 


020236142 


020235785 


36.537859 


-4.482770 


0.6206 


0.6229 


0.6218 


18.4 


425.2 


1.24 


3.9 


13 


020141586 


020141929 


36.288443 


-4.692461 


0.8310 


0.8336 


0.8323 


8.6 


425.4 


0.70 


1.6 


14 


020161356 


p020161356 


36.461850 


-4.647027 


0.9349 


0.9360 


0.9355 


10.5 


170.4 


0.43 


1.9 


15 


020162148 


020162920 


36.528063 


-4.644043 


0.6815 


0.6817 


0.6816 


8.2 


35.7 


0.03 


1.7 


16 


020170414 


020170218 


36.436325 


-4.628108 


0.3646 


0.3625 


0.3635 


14.9 


461.7 


0.46 


4.2 


17 


020182684 


020182811 


36.465562 


-4.597417 


0.7022 


0.7002 


0.7012 


14.6 


352.4 


0.24 


2.9 


18 


020198752 


020198370 


36.958101 


-4.559067 


0.9366 


0.9366 


0.9366 


12.4 


0.0 


1.22 


2.3 


19 


020281203 


020281920 


36.931009 


-4.381017 


0.9004 


0.9032 


0.9018 


11.6 


441.4 


0.99 


2.1 


20 


020383500 


020384409 


36.753843 


-4.157210 


0.5639 


0.5640 


0.5639 


16.7 


19.2 


0.64 


3.7 


21 


020323722 


020323591 


36.583688 


-4.286744 


0.9270 


0.9240 


0.9255 


12.2 


467.1 


1.21 


2.2 


22 


020336174 


p020336174 


36.883994 


-4.259516 


0.6999 


0.6988 


0.6994 


5.5 


194.1 


1.12 


1.1 


23 


020226763 


020226762 


36.771212 


-4.502420 


0.5683 


0.5662 


0.5673 


12.2 


401.7 


0.79 


2.7 


24 


020113570 


020113267 


36.553019 


-4.756721 


0.7199 


0.7218 


0.7208 


7.3 


331.0 


0.84 


1.4 


25 


020462322 


020462321 


36.730651 


-4.378359 


0.9322 


0.9323 


0.9323 


9.2 


15.5 
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Fig. 4. Postage stamps (6" x 6") and spectra of our 36 VVDS pairs selected in the bright Mg(z = 0) < -18 sam- 
ple for galaxies with r" p u " = 20h~ l kpc, Av niax = 500 km/s and AMb < 1.5. Spectra can be found either at 
http://cencosw.oamp.fr/VVDS/VVDSJ3EEP.htmli or |http://vizier.u-strasbg.fr7| 
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2.4. Accounting tor selection effects 



pair. 



To compute the total number of true pairs, we need to correct for 
three basic effects imposed by the VVDS selection function: 

1. the limiting magnitude Iab - 24 which imposes a loss of 
faint companions when we search for major mergers with 
AM B < 1.5. 

2. the spatial sampling rate and the spectroscopic success rate 
in measuring redshifts. 

3. the loss of pairs at small separations because of the ground 
based seeing limitation of the observations. 

The spectroscopic targets have been selected on the basis of 
the only magnitude criterion 17.5 < Iab ^ 24. Therefore, we 
miss companions which have an absolute magnitude fainter than 
imposed by the Iab - 24 cutoff and the AMb = 1.5 magnitude 
difference, artificially lowering the number of pairs. To take this 
into account, we compute for each galaxy a weight tJ mag (M B ,z) 
using the ratio between the comoving number densities above 
and below the magnitude cutoff (Ilbert et al., 2005). For each 
galaxy, we derive M' = M' B + AMb which corresponds to the 
maximum absolute magnitude when searching for a companion 
and M' sel (z) which corresponds to the survey limit / = 24 in the 
absolute B band at the given galaxy redshift. We then assign a 
weight for each galaxy: 



The last correction we need to apply results from the obser- 
vations which have been performed under a typical ground based 
seeing of 1 arcsecond. We correct for the increasing incomplete- 
ness to target both components of close pairs as the separation 
between them is getting smaller. Assuming a clustered distribu- 
tion of galaxies, the number of galaxy pairs should be a mono- 
tonically decreasing function of the pair separation. However, 
pairs are under-counted for separations < 2 arcseconds be- 
cause of the seeing effects. 

We derive the ratio r{6) between the observed pair count in 
the spectroscopic catalogue, N zz , over the observed pair count 
in the photometric catalogue, N pp , as a function of the angular 
separation (see Figure [5}. We apply a weight u> k g on each pair k 
using the ratio : 



" r(0 k )' 

where the mean ratio a is the probability to randomly select 
a pair, obtained at large separations. This ratio is close to 
the squared mean target sampling rate (~ 20. 2% 2 ). For large 
separations (8 > 50"), r(ff) ~ a but at small separations r(8) < a 
because of the artificial decrease of pairs due to seeing effects. 



oj' mag (M B ,z) 



1 



if AfL „ <M' , 

sup — set 



r M '~ 

%J — OQ 



<&{M)dM 






if M'„ > M' ,. 

sup set 



<b{M)dM 



We combine these weights in each pair k as co k 



■S 

J mag 



mag w mag 

, wag and bJ J mag are the weights of each galaxie in the 



pair. 

Since 25% of the field has been spectroscopicaly ob- 
served and the redshifts are not measured with 100% cer- 
tainty, we must correct for the VVDS sampling rate and red- 
shift success rate. These have been very well constrained (see 
|Ilbert et al„ 2005| ) resulting in the Target Sampling Rate (TSR) 
and the Spectroscopic Success Rate (SSR) computed as a func- 
tion of redshift. The SSR has been assumed independent of the 
galaxy type, as demonstrated to be true up to z ~ 1 in Zucca et 
al. (2006). We therefore introduce the weight (o' comp (z). For each 
galaxy, we have the information on its redshift, its apparent mag- 
nitude Iab, its spectroscopic flag and its spatial flag (whether the 
galaxy is on the field with four passes or two passes). We derive 
the completness weight as follows. 



compel 



N sel (z)\ 

gjSpectro^ 1 ' 



N: 



sel 

g, photo 



(z) ) 



• • •• 
.• • • • • 



>/« ' v v » . • ' •/ » v- • v». - 



en 

Fig. 5. Spectroscopic completeness as a function the angular pair 
separations. The line is the fit used to derive the mean correcting 
factor a which corresponds roughly to the square of the com- 
pleteness. 



The corrected number of galaxies N corr in each redshift bin 
is then: 



g 



(z) = £ 4 



where N s gspectro is the number of secure spectroscopic flag 

galaxies in the spectroscopic catalogue, and where N s vhot0 is the 

number of galaxies in the photometric catalogue. These two last 

values are estimated within the same redshift, /-band magnitude 

and N-pass area ranges based on the z, Iab and N-pass area 

values of the galaxy i. For the photometric sample, we use N c ° rr (z) — /, w p 

the photometric redshifts of Ilbert et al. (2005). Each pair k is 



comp mag' 



(2) 



The total number of pairs N"' rr is therefore computed as 



k k 

comp X ^p.mag X W fl' 



(3) 



k=\ 



Komp x ^comp where o) l comp 



therefore assigned with co pcomp 

and U) ] comp are the completeness weights of each galaxie in the 



where N g>0 bs an d Np,obs are the observed number of galaxies and 
pairs in the spectroscopic catalogue. 
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3. Evolution of the pair fraction with redshift 

3.1. Pair fraction evolution using VVDS data 

We give the total number of identified pairs as a function of the 
two separations criteria in Table 1 for the adopted AMg < 1 .5 
magnitude difference. We use equations (2) & (3) to compute 
the pair fraction f p (z) in each redshift bin as follows: 



f P (z) = 



Ajcorr 
P 



(Z) 



N c " rr (z) 



(4) 



Table[3]gives values of f p (z) for different sets of rJ" flA and Av" mx 
derived using the Mb(z) - -18 - Q(z) relation derived for the 
VVDS sample. Using the parameterization f p (z) = f P (0) x (1 + 
z) m , we fit the pair fraction measurements to compute the evolu- 
tion index m and associated poissonian errors as a function of the 
line-of-sight and projected transverse separations. These values 
are reported in Table [4] 



M B 



-18-Q(z) 



^ 



a 
o 







0.6 O.i 

redshift 

Fig. 6. Evolution of the pair fraction as a function of redshift for 
different sets of r"' ax , Av" mx = 500 km/s and galaxies brighter 
thanM B (z) = -18-g(z). 



Figure [6] shows the evolution of the pair fraction of galaxies 
with rp™* = 20, 30, 50, lOO/i" 1 kpc, A" v ,ax = 500 km/s and 
AM B < 1.5. A total of ~ 10.86 + 3.20% of galaxies with M B < 
-18 - Q(z) and r" ulx = 20h~ l kpc are in close pairs at z ~ 0.9 
compared to- 3.76+1.71% at z ~ 0.5. This leads to f p = (0.57 ± 
0.65%) x (1 + z ) 473±201 . The fraction of galaxies brighter than 
Mb — -18 - Q(z) in pairs increases significantly with redshift. 

We have investigated the dependency of the pair fraction on 
the pair separation. Increasing the separation of the two mem- 
bers of a pair both in r p and Av, the index m varies from 
4.73 + 2.01 to 2.45 + 0.11 when separations increase from 
(20/T 1 kpc, 500km/s) to (100/T 1 kpc, 2000km/s). 

Interestingly, we find a strong dependency on the limiting 
absolute magnitude of the galaxies in the pairs. Table [3] gives 
the pair fractions for different redshift, r J " ax and Av""™ using 
the M B (z = 0) < -18.77 VVDS sample and Table g] gives the 
best fit values of m and f p (z - 0). For r™ flA = 20h~ l kpc and 
Av" mx = 500 km/s, m decreases from m = 4.73 + 2.01 in the 
complete faint sample to m - 3.07 ± 1.68 for the bright sample 
with Mb{z - 0) < -18.77, implying a weaker evolution. This 
trend is seen for any separations from (20h~ l kpc, 500 km/s) 
to (100h~ l kpc, 2000 km/s). Here, we show that the pair frac- 
tion evolves faster for fainter samples. We will come back to this 
property in Section[7] 



3.2. Constraints combining low redshift pair fraction with 
VVDS estimates 

To better constrain the evolutionary parameters, the compari- 
son of high redshift data to the local value of the pair fraction 
is important. Patton et al. (2000) derived the pair fraction in a 
sample of 5426 galaxies in the SSRS2 redshift survey. Using 
close (5 < r p < 20/z kpc) dynamical (Av < 500km/s) pairs, 
they found f p (-2l < M B - 5 log h < -18) = 2.26 ± 0.52% at 
z - 0.015. We also compare our data to results from the CNOC2 
Redshift survey ( |Patton e t al., 2002) for the same magnitude se- 
lection but for a higher mean redshift : f p {—2\ < Mb - 5 log h < 
-18) = 3.21 ±0.77% at z = 0.3. 

de Propris et al (2007) derived measurements of the pair 
fraction using galaxy asymmetry and pair proximity to measure 
galaxy merger fractions for a volume limited sample of 3184 
galaxies with -21 < M B - 5 log h < -18 and 0.010 < z < 0.123 
drawn from the Millennium Galaxy Catalogue. They found a 
pair fraction of 4. 1 + 0.4% for galaxies with r p < 20h~ l kpc. 

Combining these values with our brighter sample 
(M B (z = 0) < -18 + 5log(h) ~ -18.77), we estimate 
m = 1.50 + 0.76 and f p {0) = 3.01 + 0.52. Here, we show that 
the fainter the galaxy sample is, the faster is the evolution of the 
pair fraction. 



e^ 



o 



Dh 




0.4 0.6 

redshift 

Fig. 7. Evolution of the pair fraction as a function of redshift 
adding SSRS2 (open diamond), CNOC2 (open circle) and MGC 
(open square) low redshift points to VVDS measurements (filled 
diamonds) for M B (z = 0) < -18.77. 



Figure [7] shows the best fit when combining these pair frac- 
tion measurements with our brightest sample. 

3.3. Influence of stellar mass on the pair fraction 

To identify if the evolution of the pair fraction is also dependent 
on the stellar mass of the galaxies (as a proxy for total mass), 
we applied exactly the same method as we used for the luminos- 
ity in Section [XTl but on a mass selected sample instead. Using 
masses derived in the VVDS and the evolution of the character- 
istic stellar mass, M* smr , as described in Pozzetti et al. (2007), 
we define a stellar-mass selected sample volume complete up to 
redshift ~ 1 (using an evolution parameter Qmoss(z) - -0.187xz 
to reproduce the evolution of M* lar ). Stellar masses are derived 
using a Bruzual & Chariot (2003) model and allowing bursts 
on the top of a smooth star formation history. We applied the 
same corrections described in Section 2.4 by replacing the lumi- 
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Table 3. Pair fraction (in %) for different sets of separations and redshift, with Mb < 
-18.77 - Q(z) (bright sample) using VVDS data. (See Section 3.1) 



-18 - Q(z) (faint sample) and with M B 







M7 x (z = 0) = 


-18 






I 


VI7 x (z = 0) = - 


18.77 




Av < 500km/s 




20h-'kpc 


30/r 1 kpc 


50h- l kpc 


I00hr l kpc 




20h- [ kpc 


30h- [ kpc 


50h- l kpc 


I00h- [ kpc 


z = 0.54 


3.76 ±1.71 


6.05 ±2.21 


13.78 ±3.50 


21.10 ±4.46 


z = 0.51 


2.26 ± 0.52 


7.01 ±3.00 


14.66 ± 4.54 


4.10 ±0.40 


z = 0.71 


9.43 ± 2.80 


12.52 ±3.28 


17.78 ±4.00 


40.02 ± 6.43 


z = 0.70 


3.21 ± 0.77 


12.81 ±3.65 


18.87 ±4.57 


22.03 ± 5.75 


z = 0.90 


10.86 ± 3.20 


15.30 ±3.86 


30.88 ± 5.73 


54.91 ±8.01 


z = 0.90 


4.10 ±0.40 


14.01 ± 3.70 


25.45 ± 5.20 


40.84 ± 7.25 


Av < lOOOkm/s 




20h-'kpc 


30/T 1 kpc 


50hr 1 kpc 


I00h-'kpc 




20h- l kpc 


30/T 1 kpc 


5QhT l kpc 


lOOh^kpc 


z = 0.54 


6.42 ± 2.28 


8.70 ± 2.70 


17.07 ± 3.95 


32.01 ±5.70 


z = 0.51 


7.60 ±3.14 


9.90 ± 3.64 


17.55 ± 5.04 


33.61 ±7.38 


z = 0.71 


12.44 ± 3.27 


16.27 ±3.81 


23.68 ± 4.72 


51.46 ±7.48 


z = 0.70 


10.16 ±3.20 


14.59 ±3.93 


23.14 + 5.15 


49.09 ±8.12 


z = 0.90 


13.45 ± 3.59 


19.69 ± 4.44 


38.73 ± 6.53 


66.65 ± 9.00 


z = 0.90 


11.29 ±3.28 


17.37 ±4.18 


31.29 ±5.86 


55.22 ± 8.22 


Av < 2000km/s 




20h-'kpc 


30/r 1 kpc 


50h- l kpc 


I00h-'kpc 




20/r 1 kpc 


30h- l kpc 


50h- l kpc 


I00h- l kpc 


z = 0.54 


7.22 ± 2.44 


9.50 ± 2.84 


21.34 ±4.49 


45.63 ± 7.06 


z = 0.51 


7.60 ±3.14 


9.90 ± 3.64 


20.65 ± 5.53 


41.99 ± 8.44 


z = 0.71 


12.44 ± 3.27 


17.89 ±4.02 


26.09 ± 4.99 


58.05 ± 8.05 


z = 0.70 


10.16 ±3.20 


16.45 ±4.22 


25.91 ±5.50 


55.17 ±8.73 


z = 0.90 


14.38 ±3.73 


20.62 ± 4.56 


40.40 ± 6.69 


75.65 ±9.71 


z = 0.90 


12.19 ±3.43 


18.28 ±4.30 


32.92 ± 6.04 


62.42 ± 8.86 



Table 4. Best fits parameters for m and f p (z - 0) of major mergers as a function of the dynamical parameters in the faint Mb(z 
0) < -18 sample and in the bright M B (z - 0) < -18.77 one. 





M™ a *(z = 0) = 


-18 






M™ ax (z = 0) = - 


18.77 




Av < 500km/s 


yinax 
P 


20h-'kpc 


30h-'kpc 


50h-'kpc 


lOOh-'kpc 


yinax 
P 


20h- 1 kpc 


30h- l kpc 


50h-'kpc 


lOOh-'kpc 


m 
f P (z = 0) 


4.73 ± 2.01 
0.57 ± 0.65 


4.24 ±1.36 
1.08 ± 0.83 


4.07 ± 0.95 
2.19 ± 1.18 


4.46 ±0.81 
3.27 ± 1.51 


m 
f„(.Z = 0) 


3.07 ± 1.68 
1.44 ±1.39 


3.00 ± 1.38 
2.15 ± 1.70 


2.69 ±0.16 
4.49 ± 0.42 


3.18 ±1.34 
6.22 ± 4.80 


Av < lOOOkm/s 


yinax 
P 


20/?- 'kpc 


30h- , kpc 


50h-'kpc 


lOOlr'kpc 


yinax 
P 


20h- [ kpc 


30h- l kpc 


50/7- 'kpc 


100h- l kpc 


m 
f P (z = 0) 


3.37 ± 1.52 
1.66 ± 1.43 


3.77 ± 1.13 
1.86 ± 1.18 


4.05 ± 0.55 
2.84 ± 0.88 


3.46 ± 0.54 
7.46 ± 2.26 


m 
f„(.z = 0) 


1.81 ±0.51 
3.62 ± 1.04 


2.57 ± 0.57 
3.41 ± 1.11 


2.80 ±0.10 
5.15 ±0.31 


2.25 ± 0.70 
13.4 ±5.36 


Av < 2000km/s 


yinax 
P 


20h-'kpc 


30h-'kpc 


50h-'kpc 


I00h-'kpc 


jnax 
P 


20/r 1 kpc 


30h- l kpc 


50h-'kpc 


1 00/r l kpc 


m 
f P (Z = 0) 


3.19 ± 1.04 
1.96 ± 1.15 


3.56 ±1.27 
2.24 ± 1.61 


3.18 ±0.75 
5.10 ±2.14 


2.45 ±0.11 
15.7 ±0.98 


m 
f P (z = 0) 


2.21 ±0.29 
2.98 ± 0.49 


2.71 ± 1.08 
3.36 ± 2.08 


2.25 ± 0.06 
7.72 ± 0.27 


1.83 ± 0.40 
19.6 ±4.43 



nosity function by the mass function. We define a major pair 
via the ratio M1/M2 of stellar masses, and select pairs with 
Mi IM2 < 4 corresponding roughly to a luminosity selected sam- 
ple with AMb < 1.5 mag. 

We divided our sample in different sub-samples: one with 
log(M/M e ) > 9.5 (106 pairs), one with log(M/M e ) > 10 (77 
pairs) and one with log(M/M G ) > 10.5 (37 pairs) with sep- 
arations Av < 500 km/s and increasing the projected separa- 
tion tor[" ax - 100/tT 1 kpc for better statistics, as shown in 
Figure [8f Figure [9] shows the evolution of the pair fraction in 
those different mass sub-samples. For low mass galaxies with 
log(M/M e ) > 9.5, m = 3.13 + 1.54 and f p (z = 0) = 3.90 + 3.42. 
For intermediate mass galaxies with log{M/M Q ) > 10, m - 
2.04+1.65 and f p (z — 0) = 7.28±6.81.Formassive galaxies with 
log(M/M Q ) > 10.5, m = 0.52 + 2.07 and f p (z = 0) = 16.7 ±19.5. 
We see a flatter evolution as we select more massive galaxies. 
It is therefore apparent that intermediate or low mass galaxies 



are responsible for most of the evolution of the pair fraction and 
merger rate. 

4. Physical properties of galaxy pairs 

4. 1 . Spectro-photometric properties 

One of the expected effect of a merging or close interaction of 
galaxies is an increase in the star formation rate of the system. 
We evaluate here if our sample of pairs has a stronger star for- 
mation rate than the global population by studying the rest-frame 
[OH] 3727 A equivalent widths (EW) as a function of projected 
separation within a given Av" rav (500 km/s). £W[0/7] were 
derived using the plate fit software (Lamarei lle et al., 2 006), 
which performs a continuum fit to the observed spectra using 
template fitting. It enables an unbiased measurement of the in- 
tensities of absorption and emission lines. For each pair, we pro- 
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Table 5. Spectral types of pairs for the less and the most massive selected samples. Fractions are given in brackets. 
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N D , 



Early-type pairs Late-type pairs Mixed type pairs 



log(M/M ) > 9.5 106 31(29.2%) 53(50.0%) 22(20.8%) 

log(M/M e ) > 10 77 29 (37.7%) 29 (37.7%) 19 (24.6%) 

log(M/M e ) > 10.5 37 22 (59.5%) 6 (16.2%) 9 (24.3%) 
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Fig. 8. The three sub-samples defined to study the influence of 
the mass on the pair fraction using pairs with r J " ax = 100/T 1 kpc. 
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Fig. 9. Evolution of the pair fraction for different sub-samples 
with different stellar mass limits log(M/M Q ) > 9.5 (circles), 
log(M/M Q ) > 10 (squares) and log(M/M Q ) > 10.5 (empty di- 
amonds). 



duced the mean EW[OII] by summing the individual EW[OII] 
of each galaxy, assuming EW[OII] - if the line is not de- 
tected. Using only galaxies for which the [OH] line has been 
detected, the mean EW[OII] is larger, on average, for the clos- 
est pairs with EW[OII] = 46.7 + 4.35 for r™ flJC < 20/T 1 kpc, 
EW[OII] = 40.5 + 3.78 for r" flA < 50/r 1 kpc, and EW[OII] = 
36.5 ± 3.12 for r" flA < 100/T 1 kpc indicating a 25.9 ± 4.10% 
increase in EW[OII] at small separations (see Figure ITOb . We 
perform the same estimation using also galaxies for which the 
[Oil] line has not been detected (EW[OII] - 0). Both samples 
show an increase of the mean EW[OII] at small projected sep- 
arations, extending to higher redshifts results of Woods et al. 
(2006) in the local CFA2 sample. We conclude that star forma- 



tion is enhanced in close merging systems at the mean redshift 
< z >= 0.76 of our sample. 
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Fig. 10. Mean EW[OH] for pairs with Av' 



500 km/s as a 



function of rT*. We present results using only galaxies where an 
[OH] line has been detected (filled circles) and using all galax- 
ies, including those where the line has not been detected (empty 
circles). Thin lines represent the mean value of EW[OII] for the 
two sub-samples fitted on 80 < r™" < 150h~ l kpc. 



4.2. Spectral types of galaxies in pairs 

In this Section, we compare the spectral properties of galax- 
ies in dynamical pairs with field galaxies. For each galaxy 
in the VVDS, the spectral type has been derived on the ba- 
sis of the template fitting of the rest-frame multi-/l photometry 
( |Zucca et~ al., 2006). Galaxies were classified in type 1 (E/S0), 
type 2 (Early spiral), type 3 (Late spiral) and type 4 (Irregular). 
Therefore, for each pair, we know the spectral types of both the 
primary galaxy and its companion(s). 

We have investigated which galaxy types are involved in a 
pair as cosmic time evolves. We classified each pair with a flag 
(X - X) where X is the spectral type of each pair member. For in- 
stance 'dry mergers' with the merging of two early-type galaxies 
are classified as type (1 - 1). We consider all the permutations 
between these four types. We classify as 'early-type' pairs, pairs 
with flags (1 - 1), (1 -2) and (2 - 2), late-type pairs the pairs 
with flags (3 - 3), (3 - 4), and (4 - 4), and mixed type pairs those 
with flags (1 - 3), (1 - 4), (2 - 3), and (2 - 4). Table gives 
the fraction of these different classes in the r"' ax = 100/i~' kpc 
pair sample. The fraction of pairs involving only E/SO galaxies 
increases from 3.0% at z ~ 0.9 to 11.8% at z ~ 0.5, the frac- 
tion of pairs involving at least one E/SO increases from 22.4% 
at z ~ 0.9 to 29.4% at z ~ 0.5, while the vast majority of pairs in- 
volving at least one late-spiral or Irr galaxy represents a fraction 
decreasing from 83.4% at z ~ 0.9 to 76.5% at z ~ 0.5. 
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Table 6. Fraction of pairs vs. the spectral classes of each galaxy in the pair for redshift z ~ 0.5 and z ~ 0.9. Pairs with Mb < 
-18 - Q(z) and r J " ax = lOOkpc/h are considered (202 pairs in total). 



Classification Fraction at z ~ 0.5 Fraction at z ~ 0.9 



2E/SO 

1 E/SO involved 

1 E/SO or 1 early-Sp involved 



(1-1) 
(1-X) 

(l-X) + (2-X) 



11.8% 
29.4% 
47.1% 



2Irr 

1 Irr involved 

1 Irr or 1 late-Sp involved 



(4-4) 

(4-X) 

(3 - X) + (4 - X) 



17.6% 
47.1% 
76.5% 



3.0% 
22.4% 
46.3% 



16.4% 
47.8% 
83.4% 



;r 



a 
o 



£ 




0.4 0.6 

redshift 

Fig. 11. Fraction of early-type and late-type galaxies in the un- 
derlying selected sample brighter than Mb — -18 - Q(z) (filled 
squares) and in the underlying selected sample brighter than 
Mb - -18.77 - Q(z) (empty squares) as a function of redshift. 
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Fig. 12. Fraction of early- (circles), mixed- (diamonds) and late- 
type (squares) pairs in the selected sample brighter than Mb - 
-18 - Q(z) as a function of redshift. 



Figure QT| shows the evolution of the fraction of early-type, 
and late-type galaxies in two different samples : one brighter 
than Mb = -18 - Q(z) (faint sample) and one brighter than 
Mb = -18.77 - Q(z) (bright sample). In the faint sample, the 
population is dominated by late-type galaxies at all redshifts. In 
the bright sample, late-type galaxies dominate between z ~ 0.4 
and z ~ 1, and early-type galaxies become dominant between 
z ~ 0.1 and z ~ 0.3. Early-type galaxies represent only one 



third of the sample at z ~ 1, but about two third at z ~ 0.1. 
Figure [12] shows the fraction of early, mixed and late-type pairs 
with M B < -18 - Q(z) as a function of redshift. At z ~ 0.9, 15% 
(55%) of these pairs are early(late) type pairs whereas at z ~ 0.5, 
25% (50%) of these pairs are early (late) type pairs following the 
same trend as the underlying sample of galaxies. 

FigureQ~3]shows the pair fraction as a function of redshift, for 
3 classes of pairs type compared to the global pair fraction: two 
early-type galaxies, one early and one late component, two late- 
types, using a sample with r"' ax - 100/i~' kpc. The early-type 
pair fraction evolves slowly with redshift with m = 1.44 + 0.93. 
On the contrary, the late-early and late-type pair fractions evolve 
strongly with redshift, with to = 5. 16 ±2.56 and to = 4.74 ±0.81 
respectively. 

Table [5] gives the distributions of pairs as a function of 
the stellar mass selection and spectral types of the pairs. The 
log(M/M Q ) > 9.5 sub-sample is dominated by late-type pairs 
(50%) while the log(M/M e ) > 10.5 sub-sample is dominated 
by early-type pairs (59.5%) (see Section 4.2). We conclude that 
most of the pair fraction evolution is coming from lower mass 
late-type or mixed-type pairs. 
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Fig. 13. Evolution of the global pair fraction in the faint sample 
(M B (z = 0) < -18) with r™ flA = lOO/i" 1 kpc (pink) and con- 
tribution of early-types (red), mixed types (cyan) and late-type 
(blue). 



5. Evolution of the merger rate 

Knowing the pair fraction, we derive the merger rate i.e. the 
number of mergers per unit time and per comoving volume. This 
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rate can be expressed as 

(\rcorr _ jycorr \ 
V P triplets } 



N mg (z) = C, 



ni: 



x n(z) x T. 



"•!■! 



where C mg stands for the fraction of galaxies in close pairs that 
will undergo a merger within the time T mg and n(z) is the co- 
moving number density of galaxies. The best way to estimate 
these values is to use simulations to follow the merging his- 
tory of galaxies with different masses. We take results from the 
Millennium simulations (Kitzbichler & White, 2008) to estimate 
the merging time-scale T mg (r"^ x ,z), written as follows: 



— 1/2 



T (r"' ax r U2 + fi(rT x ) x z + fiiC*) * (bgM. - 10) 



We computed To, f\ and fi for r^ ax — 

(20, 30, 50 and lOO)^ 1 kpc in the case of Av max = 500 km/s. 
Following Lin et al. (2008), the probability for a pair to merge 
in the given time-scale T mg is assumed constant, C mg - 0.6, 
independent of the separation r™ ax . As a proxy for total mass, 
we use the evolution of the characteristic stellar mass M* tars as 
derived in the VVDS flPozzetti et al., 2007| l. Figure [TU shows 
the change in the galaxy merging time-scale with redshift and 



merger rate increases from ~ 12.3x10 4 
ers h? Mpc~ 3 Gyr 1 from z = 0.5 to z 



to- 19.4xl0~ 4 merg- 
= 0.9. The merger rate 



evolves as N, 



(5) and N mg (z 



mg\Z) — -'»/ 

= 0) = (4.96 



mg (z = 0)x(l+z) m »« withm mg = 2.20+0.77 



2.07) x 10~ 4 . Table lists the values 
(z = 0) for different sets of sepa- 
rations and we plot in Figure [L5]the evolution of the merger rate 
30, 50, 100/T 1 kpc and Av max = 500km/s. 



of the parameters m mg and N mg yi, 



for r" p mx = 20 



5«. 



a, 






B 

C3 
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Fig. 15. Evolution of the Mb < -18 - Q(z) galaxy merger rate 
for different sets of t J ™ x with AV" ax - 500km/ s. 
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Fig. 14. Evolution with redshift of the merging time-scale in Gyr 



as a function of r™ 



The time-scales are found to be higher than the standard 
assumption that half of the pairs with r^ ax = 20h~ l kpc 
undergo a merger in half a Giga-year ( |Patton et al., 2 000. 
IPatton et al., 20021 |Lin et al., 2004) 1. Using the Kitzblicher and 
White (2008) prescription, we find that even for the closest pairs 
the merging time-scale is 1.5 times higher than assumptions pre- 
viously used in the literature. 

The merger rate should be an "absolute value", independent 
of r" rav and Av since we take into account the merging time- 
scales corresponding to different pair separations. To check that 
the merger rate does not depend on the adopted value of r™ ax 
and Av, we have computed the merger rate for different sets of 



Using this merger rate evolution parametrisation, we esti- 
mate the fraction of present day galaxies f rem , that have under- 
gone one major merger (Patton et al., 2002) since z ~ 1. 



f n 



N , 

-'-ni^ 



fm g (Zj) 



SfmgiZj) 



(6) 



where Zj corresponds to a lookback time of t — j x T mg and f mg is 
the fraction of galaxies that undergo a merger. We use the merger 
rate evolution derived with r™ ax = 20/i~' kpc and lookback times 
derived using 



^lookback — TT 



dz' 



H oJo (l+z')V" m (l+z') 3 + ^ 



(7) 



based on a mean merging time-scale of 0.75 Gyr (corresponding 
to r™ flA = 20h~ l kpc). We find that 8% of present day galaxies 
brighter than Mb = -18 - Q(z) have undergone a major merger 
since z ~ 0.4, , while 22% have done so since z ~ 0.9. 

We have also computed the merger rate for two different lu- 
minosities using pairs with r"' flA = \00h~ l kpc. A similar trend to 
the pair fraction is observed: for galaxies with Mb(z - 0) < -18, 
we find m. 



a IK 



2.20 + 0.77, while for brighter galaxies with 
M B (z - 0) < -18.77 we find m mg = 1.60 + 1.83 using only 
VVDS data. For the same limiting magnitude and using the 
merger rate measured by de Propris et al. (2007) to constrain 
the low redshift end, m mg = 1.57 + 0.44. 

Similarly, we have computed the merger rate for different 
mass selected samples as defined in Section 3.3 using r p wx = 

r"' ax with Av < 500km/ s; results are presented in Table and I00h l kpc. For the less massive sample {log(M/M ) > 9.5) 



plotted in Figure [15] The merger rate values are in good agree- 
ment, both in slope and normalization, for different sets of pro- 
jected separations. This is a good indication of the robustness of 
the method. In the following, we use values of the merger rate 
with r"' ax - 100h~ l kpc for better statistics, when necessary. The 



= 2.38 + 1.57 with N mg (z = 0) = (3.56 + 3.17) x 
10~ 4 mergers h 3 Mpc~ 3 Gyr~ l , while for the intermediate sam- 
ple (log{M/M G ) > 10), m mg = 1.27 + 1.67 with N mg (z = 
0) = (2.75 + 2.61) x 10~ 4 mergers h 3 Mpc- 3 Gyr\ as shown 
in Figure [T71 
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Table 7. Merger rate values for different sets of parameter and redshift in units of 10 4 mergers Mpc ^Gyr ' for the bright sample 
(M" B mx {z = 0) = -18) for galaxies with AM B < 1.5. 



Av < 500 km/s 



20h- [ kpc 



30/r 1 kpc 



50hr l kpc 



I00h- [ kpc 



Zmean 
Zmean 
<-mean 



: 0.5124 8.17 ±3.87 
: 0.6952 15.60 ± 5.97 
: 0.8989 14.08 ± 6.47 



8.95 ± 3.39 
14.10 ±4.73 
13.48 ± 5.27 



14.43 ± 3.56 
14.18 ±3.85 
19.31 ±5.21 



12.32 ± 2.41 
17.91 ±3.29 
19.37 ± 3.90 



Av < 1000 km/s 



20h- [ kpc 



30h- [ kpc 



50h~ l kpc 



I00h- [ kpc 



Zmean 
Zmean 
■(•mean 



0.5124 13.94 ±5.15 

0.6952 20.59 ± 6.95 

; 0.8989 17.44 ± 7.26 



12.89 ±4.13 
18.32 ±5.47 
17.35 ± 6.05 



17.92 ±4.01 

18.93 ±4.53 
24.27 ± 5.91 



18.69 ± 3.07 
23.02 ±3.81 
23.51 ±4.37 



Av < 2000 km/s 



20h- l kpc 



30h-'kpc 



50h-'kpc lOOh^kpc 



Zmean 
Zmean 
■^mean 



0.5124 15.67 ±5.49 
0.6952 20.59 ± 6.95 
0.8989 18.65 ±7.52 



14.07 ±4.33 
20.14 ±5.77 
18.17 ±6.20 



22.58 ±4.55 
21.03 ±4.78 
25.53 ±6.05 



27.17 ±3.78 
26.48 ±4.10 
27.21 ±4.70 



Table 8. Best fit parameters of m mg and N mg (z - 0) for major mergers as a function of the dynamical parameters for the faint sample 
and for galaxies with AM# < 1.5. . 



Av < 500 km/s 



N mg (z = 0) x 10" 



lOh^kpc 

2.63 ± 1.96 
2.93 ±3.14 



30hr l kpc 
2.01 ± 1.32 
4.03 ± 2.88 



50ft" 1 kpc 

1.33 ±0.98 
7.70 ± 4.05 



Av < 1000 km/s 



20h- [ kpc 30h' l kpc 50h- [ kpc 



N„,Jz = 0) x 10" 



Av < 2000 km/s 



N m Jz = 0) x 10" 



0.93 ± 1.00 
11.0 ±5.89 



1.31 ±1.23 
8.49 ± 5.61 



0.49 ± 0.78 
17.5 ± 7.24 



lOOh^kpc 

2.20 ± 0.77 
4.96 ± 2.07 



I00h- [ kpc 



1.19 ±1.48 1.50 ±1.08 1.42 ±0.58 1.15 ±0.50 
8.98 ±7.12 7.14 ±4.16 9.40 ± 2.94 11.7 ±3.10 



20h- l kpc 30h~ l kpc 50h- l kpc I00h~ l kpc 



-0.02 ±0.15 

27.2 ±2.13 



We see a change in the evolution of the merger rate as we go 
to the highest masses. First, the number of less massive merging 
events (log(M/M Q ) > 9.5) is greater than the number of high 
mass merging events (log{M/M e ) > 10.5). Then we see a flat- 
tening of the evolution of the merger rate as we go to higher mass 
galaxies, confirming that the evolution of the major merger rate 
is mainly due to the less massive galaxy population. 

6. Stellar mass involved in mergers 

We estimate the fraction of the total stellar mass involved in a 
merger process, f M > (z), since z ~ 1 as a function of redshift as 



/m'(z) 



Ml 



(z) x N m Jz) x T bin (z) 



M*(z) x n(z) 



(8) 



where n(z) is the comoving number density of galaxies, N mg (z) 
is the number of mergers per unit of time and per comoving vol- 
ume, M* (z) - N ' + f J is the mean stellar mass involved in 
a merger process, Tbin(z) is the elapsed time corresponding to 



the considered redshift bin and M* ot (z) is the total stellar mass 
in the redshift interval. To extrapolate the values of the stellar 
mass densities at z ~ 0.1, we assumed a constant stellar mass 
density below z - 0.4. This assumption is consistent with the 
evolution of p„ reported in Pozzetti et al. (2007). We show in 
Figure[18]that around 25% of the stellar mass contained in galax- 
ies with log(M/M G ) > 9.5 at z ~ 0.1 have experienced a merger 
since z ~ 1 while this fraction is about 20% for galaxies with 
log{M/M Q ) > 10. One can identify two trends: the fraction of 
the stellar mass density coming from the merging process shows 
a rise of about 24% from z ~ 0.9 down to z ~ 0.1 for the less 
massive population, whereas it stays roughly constant at about 
20% for the most massive galaxies. 



7. Summary and Discussion 

Our results can be summarized as follows: 

(i) We find that 3.8 + 1.7, 9.4±2.8,and 10.9+3.2% of galaxies 
with M B {z) < -18 - Q(z) at z ~ 0.5, 0.7 and 0.9 respectively, are 
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Fig. 16. Evolution of the merger rate for different luminosity 
ranges. For the brightest sample (Mb{z — 0) < -18.77), we 
add results from de Propris et al. (2007) (empty square) 
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Fig. 17. Evolution of the merger rate for different mass ranges 
using r™ flA = lOOh^kpc. From top to bottom : log(M/M Q ) > 9.5 
(diamonds), log(M/M Q ) > 10 (circles) and log(M/M Q ) > 10.5 
(squares) 
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Fig. 18. The fraction of stellar mass density involved in a merger 
process since z ~ 1 as a function of redshift for different mass 
selected sub-samples. 



in pairs of galaxies with luminosities AM# < 1.5 and separations 
less than 20/i~' kpc. 



(ii) The evolution of the pair fraction with redshift is strongly 
dependent on the absolute luminosity or stellar mass of the 
brighter galaxy in the pair: it evolves more slowly for brighter or 
more massive galaxies than for faint galaxies. Using the VVDS 
alone, the pair fraction of galaxies with Mb(z) < -18 - Q{z) 
is found to strongly evolves with redshift as oc (1 + f)" 1 with 
m - 4.46 ± 0.81 for separations of (100/T 1 kpc, 500km/ s), 
while for brighter galaxies with Mb(z) < -18.77 - Q(z), we 
find a slower evolution with m = 3.18 + 1.34. Combining 
VVDS data with low redshift measurements from de Propris et 
al. (2007), Patton et al. (2000) and Patton et al. (2002), and tak- 
ing r J ™ x - 20/i~' kpc, we similarly find m = 1.50 + 0.76 for 

bright galaxies with M B (z = 0) < -18 + 5log(h) 18.77 and 

m - 4.73+2.01 for the fainter M B (z - 0) < -18 sample. In addi- 
tion, the evolution of the pair fraction is found to be stronger with 
m = 3. 13 + 1.54 for less massive galaxies with log{M/M e ) > 9.5, 
than for more massive galaxies with log{M/M Q ) > 10 for which 
we find m - 2.04 ± 1 .65. Low mass pairs are therefore contribut- 
ing more to the evolution of the pair fraction than high mass 
pairs. 

(iii) The star formation rate of close pairs is enhanced at sep- 
arations r p < 150/i _1 kpc. We find that the mean EW(OII) in 
close pairs are larger by 26 ± 4% than the one derived for galax- 
ies with larger separations. 

(iv) The evolution of the pair fraction is stronger for late- 
type pairs with mi ale - 4.74 ±0.81, than for early-type pairs with 
m-eariy = 1-44 ± 0.93. Late-type pairs are therefore contributing 
significantly more to the observed evolution of the pair fraction 
than early-type pairs in our Iab S 24 sample. 

(v) Using the merging timescale from Kitzbichler & White 
(2008), we find that the merger rate increases from ~ 12.3 x 
10~ 4 to ~ 19.4 x 10~ 4 mergers h 3 Mpc^ Gyr 1 from z - 0.5 
to z — 0.9. The merger rate of galaxies with Mb(z) < -18 - 
Q(z) evolves as N mg = (4.96 ± 2.07) x 10~ 4 ) x (1 + z ) 2 - 20±() - 77 . 
Similarly to the pair fraction, we find that the merger rate evolves 
faster for fainter or less massive galaxies, with m mg - 2.20 + 0.77 
and 2.38 + 1.57 respectively, than for brighter or more massive 
galaxies with m, 



nix 



1.57 + 0.44 and 1.27 + 1.67 respectively. 
The merger rate is evolving more strongly for late-type mergers 
than for early-type mergers. 

We conclude that the observed evolution of the pair fraction 
and merger rate in our Iab ^ 24 sample is mostly driven by 
low mass late-type galaxies, while the pair fraction and merger 
rate of high mass early-type galaxies remains roughly constant 
since z ~ 1 . Therefore, the pair fraction or the merger rate are 
not universal numbers but rather are dependent on the luminos- 
ity or stellar mass, and on the spectral type of galaxies involved. 
Our finding that bright or massive galaxies experience a lower 
merger rate and a lower evolution of the merger rate extends 
to higher redshifts the results found in the local Universe by 
Patton & Atfield (2008). Taking into account this pair fraction 
and merger rate dependancy on galaxy luminosity and spectral 
type offers a first step to reconcile apparently inconsistent ob- 
servations. Lotz et al. (2008) find a slow or no evolution of the 
merger rate and claim that they disagree with previous studies. 
When taking into account that their result is derived from bright 
Mb < -19.94 - 1.3 X z galaxies, their result is consistent with 
other studies like Conselice et al. (2003) or Le Fevre et al. (2000) 
which have been analysing fainter samples. 

The dependency of the merger rate and its evolution on lu- 
minosity or stellar mass is indeed a prediction from the latest 
simulations using advanced semi-analytic models as described 
in Kitzbichler and White (2008). At the limiting magnitudes or 
stellar masses of our sample, Kitzbichler and White (2008) pre- 
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diet that the merger rate decreases and evolves more slowly for 
galaxy samples with increasing luminosity or stellar mass, simi- 
lar to the trend observed in our sample. 

The star formation rate is significantly enhanced in merging 
pairs with a net star formation increase of ~ 25% for these galax- 
ies. Nevertheless, it accounts for only 12% to 3% of the global 
galaxy population from redshift z ~ 1 to z ~ which is not suffi- 
cient to counteract the strong fading of the global star formation 
rate observed since z ~ 1 . This may indicate that the gas reser- 
voir of massive and intermediate mass galaxies has already been 
depleted at redshifts z ~ 1, in agreement with their observed 
peak in star formation at z ~ 3.5 (e.g. Tresse et al., 2007). It is 
then apparent that the decreasing SFR since z ~ 1 is regulated by 
other physical processes like gas availability in the intergalactic 
medium, or feedback. 

Major merging events are largely dominated by pairs of late 
or mixed type galaxies, but while early-type mergers represent 
about 15% of the merging events of bright galaxies at z ~ 1, 
they become approximately 25% of all mergers at z ~ 0.5, 
which is in good agreement with previous results on dry mergers 
(e.g. |Lin et al., 2008| . This indicates that major mergers are effi- 
cient in lowering the number density of intermediate mass late- 
type galaxies to build up more early-type galaxies. We confirm 
that merging is one of the important physical processes driving 
galaxy evolution, with the observed galaxy merger rate undoubt- 
edly closely linked to the hierarchical build up of dark matter 
galaxy halos, with a rapid mass accretion phase of massive halos 
since z ~ 1 (Abbas et al., 2008). Our finding that ~ 20% of the 
stellar mass in present day massive galaxies has experienced a 
major merger since z ~ 1 is an indication that major mergers are 
significantly contributing to the observed evolution of the stel- 
lar mass density since z ~ 1 (Bundy et al., 2005; Arnouts et al., 
2007). 
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